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Abstract 



We study the coupling of Abelian gauge theories to four- dimensional simplicial quan- 
tum gravity. The gauge fields live on dual links. This is the correct formulation if we 
want to compare the effect of gauge fields on geometry with similar effects studied so 
far for scalar fields. It shows that gauge fields couple equally weakly to geometry as 
scalar fields, and it offers an understanding of the relation between measure factors 
and Abelian gauge fields observed so-far. 
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1 Introduction 



A path integral formulation of quantum gravity includes in its simplest version 
an integration over all four-dimensional geometries of a given fixed topology, the 
weight of each geometry being the Boltzmann weight of a suitable action. An 
attempt to implement such a path integral representation is known as dynamical 
triangulation. The integration over geometries is approximated by a summation 
over all triangulations constructed from equilateral simplices of side- length a, and 
the continuum limit is obtained by letting a — > for a suitable choice of bare coupling 
constants in the action |], |J. 

This model has been analyzed in a series of papers, both with and without cou- 
pling to scalar fields. Originally the hope was that an observed second order phase 
transition could be used to define a non-perturbative theory of quantum gravity, 
but a closer examination revealed that the transition was (weakly) first order ||. 
Motivated by an effective theory of quantum gravity, which showed that the gauge 
fields coupled very strongly to the infrared sector of gravity (3J , the coupling between 
Abelian gauge fields and simplicial quantum gravity was introduced and studied us- 
ing the non-compact version of the gauge field IJ. We shall refer to this version of 
the model as AGM (Abelian gauge model). Indeed, for the first time a significant 
coupling between matter and gravity was observed, and a new phase, named the 
"crinkled" phase was observed. It appeared to be different from the so-called "crum- 
pled" and "elongated" phases observed that far, and which both seemed irrelevant 
for a continuum limit of quantum gravity. The crinkled phase had a negative en- 
tropy exponent 7 and a fractal dimension du ~ 4. However, it was later shown that 
most of the properties of this new phase could be obtained by a simple change in the 
measure of the path integral || , and since such a change is of ultra-local nature it is 
unlikely to provide a faithful representation of the ideas relating the trace anomaly 
to the infrared behaviour of quantum gravity, assuming conformal invariance. 

While the change in measure seems to capture most of the gravitational physics 
associated with the coupling of Abelian matter fields to simplicial quantum gravity, 
the two theories are of course not equivalent, as emphasized in [^, [?J . One can explore 
the full consequences of the Abelian gauge fields coupled to simplicial quantum 
gravity by Monte Carlo simulations. Such studies were performed in |J and a few 
ambiguities as to the nature of the transition to the crinkled phase and of the crinkled 
phase itself remained. For this reason and as a matter of principle it is valuable to 
have an independent simulation of the system. 

In the non-compact formulation the Abelian gauge action is Gaussian. There 
exist well-known methods relating low- dimensional Gaussian theories on a direct 
lattice to their dual versions. The methods developed on regular lattices work equally 
well on dynamical lattices and we shall make use of the concept of duality to relate 
theories defined on direct and dual lattices. This relation will highlight why the 
originally chosen way of coupling gauge fields to geometry produced results almost 
identical to a simple measure term. It will show that this is the result of a somewhat 
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unfortunate choice of action, and that the "real" Abelian degrees of freedom couple 
to simplicial quantum gravity more or less as ordinary massless scalar fields, i.e. they 
do not influence the phase structure of pure simplicial quantum gravity, at least for 
the number of independent Abelian fields used so far. From this point of view 
the crinkled phase observed so far is entirely due to the insertion of an ultra-local 
measure factor. Further, the way we implement the coupling between the gauge 
fields and geometry allows us to introduce new effective actions which might change 
the phase diagram of matter coupled to gravity. 

Another possibility, not studied in this paper would be to go to the compact 
formulation of the gauge action. In this case, at least in principle, one may expect 
a phase transition between the confined and deconfined phases. Numerical results 
for the three-dimensional case suggest however that such a transition disappears, at 
least for a single U(l) field ||. Whether some new physics may emerge for more 
fields or for the non- Abelian gauge fields remains an open problem. 

2 The model 

Let T be an abstract triangulation of a four-dimensional manifold. Let T* be the sim- 
plicial complex dual to T and denote the number of (i-dimensional sub-simplices 
in T. N 4 is the total volume of the manifold. We shall compare two possible methods 
of introducing the (non compact) Abelian gauge fields on such a manifold. 

The first possibility, which we call the T implementation, would be to put 
the (complex) matter fields at the vertices in T and the U(l) gauge potentials 
Au = —Aji on the (oriented) links (IJ) of the triangulation T, with /, J number- 
ing the vertices of the triangulation. At each vertex I we can perform the gauge 
transformation \i an d the gauge fields transforms as 

A u -> A u + X i- Xj- (1) 

Let us now choose a spanning tree in the complex T, i.e. a connected graph of links 
with no loops. The number of links in the spanning tree is Nq — 1. Making use 
of the gauge transformation ([!]) we can perform the gauge fixing of all the gauge 
potentials on the spanning tree, reducing them to zero. The resulting number of 
degrees of freedom is N% — No + 1. In principle there is nothing wrong with this 
implementation where the matter fields live on the vertices. However, the number 
of the degrees of freedom both of the matter fields and the gauge fields has no 
simple relation to the volume of the manifold. For the scalar fields this number is 
proportional to the number of vertices Nq (and not N4), which can grow as N£, 
where a < 1. In cases where one has a sensible thermodynamic limit it should 
not matter which way the coupling between matter and gravity is done, but since 
"pathological" triangulations can occur, and even dominate for some choices of 
gravitational coupling constants, this choice of implementation may have troubles 
with defining the extensive quantities, like the average action etc. 
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The alternative possibility, which we call the T* implementation, is to choose 
the £7(1) gauge potentials A*j to live on the oriented links of the dual graphs, 
connecting simplices with labels i and j. Each dual link corresponds to an interface 
between the neighbouring simplices. As usual A* { = —A*j and there are N 3 = -A^ 
such fields. This formulation is in accordance with most simulations of matter 
fields coupled to gravity. The matter fields are usually associated with the four- 
simplices (are located in the "centers" of the four-simplices) and their derivatives 
are calculated as differences formed between the values of the fields in neighbouring 
four-simplices, i.e. they are associated with the dual links. The reason for choosing 
such an implementation, apart from being convenient in numerical simulations, is 
that the naive counting of field degrees of freedom will be proportional to the four- 
volume A/4. This property holds also for the gauge degrees of freedom. At each 
simplex we can perform a gauge transformation \l and the gauge fields transform 
as 

^-^ + Xj-Xj- (2) 

Let us now choose a spanning tree in the complex T* . The number of links in the 
spanning tree is A/4 — 1. One can gauge A*, = for all links in the spanning tree. 
This leaves us with A3 — N4 + 1 = |AT 4 + 1 physical degrees of freedom for the gauge 
field. The nice feature is that this number depends only on the volume A/4. 

In both implementations the independent degrees of freedom can be parametrized 
by the gauge invariant (oriented) plaquette variables. In the T implementation these 
are oriented triangles tjjK and the plaquette variables satisfy: 

P(t JJK ) = A u + A JK + A KI . (3) 

The plaquettes are not independent variables. Because of the Abelian nature of 
the gauge fields Pitux) are unchanged under the cyclic permutation of J, J, K and 
change sign when the orientation of the plaquette is inversed. For each triangula- 
tion one has to define the positive orientation for the plaquette variables to avoid 
the double counting. By construction the flux through any closed two-dimensional 
surface is zero. The number of independent such surfaces is (in case the topology is 
that of the 4-sphere, which we will assume) 

N(S) =N 3 -N A + 1. (4) 

This is the number of independent Bianchi identities we can write down. By Euler's 
relation for the 4-sphere it leaves precisely A" 4 — N + 1 independent plaquettes, and 
the Jacobian for changing from independent gauge fields to plaquette variables is 
one [0: 

/ n dA i n *(4o F ( p (A)) = I n *p® n *e pm nn (5) 

J leT I'eST t s es tes 

In eq. @ ST denotes a spanning tree in T and S denotes the set of N(S) independent 
two-dimensional surfaces where the Bianchi identities are enforced. I and t are 
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respectively links and triangles of the manifold and we assume that the definition 
of the positive orientation of links and triangles was uniquely chosen. Notice that 
in this implementation each elementary gauge loop (plaquette) has a length three, 
however each gauge potential contributes to a number of neighbouring triangles 
which is a non-trivial local geometric characteristic of the system. 

Similar discussion can be made in the T* implementation. The dual plaquettes 
constructed from a .D- dimensional triangulation can be labelled by the (D — 2)- 
dimensional sub-simplices in the D-dimensional triangulation, which are "encircled" 
by the associated dual plaquette. In four dimensions we shall label dual plaquettes 
by the triangles they encircle. The geometric properties of the two gauge-invariant 
objects P(U) and P*(ti) are however different: The length of the plaquette P*(ti) is 
now o(ti), where o(ti) is the order of the triangle t iy or the number of four-dimensional 
simplices sharing this triangle. As before we have to choose a positive orientation, 
and the plaquette variable is unchanged when the corresponding loop is rotated and 
changes sign, when it is inverted. Each gauge potential A*- contributes to exactly 
four dual plaquette variables. Like in the T implementation we can parametrize the 
system by the plaquette variables P*(ti), satisfying a set of the Bianchi identities. 

To this point we did not specify the model we want to consider. Let us start 
with the T* implementation. The plaquette variables are the sums of gauge fields 
associated with the oriented (dual) links 'encircling' the triangles tt in T: 

^)=£^*(a (6) 

a 

where the orientation of the (dual) links follows that of the plaquette and as discussed 
above U £ T can be viewed as dual to the plaquette. To keep the discussion general 
let us at this point assume only that the gauge action is Gaussian and assumes a 
form: 

S* g (T,A*) = J2p i P* 2 (t t ), (7) 

i 

where fa are positive, local, geometry-dependent coefficients. We shall argue below 
what seems to be the most natural choice of fa. 

With this notation the partition function can be written as 

Z = ^-±-e-< N *W +l * N *W /lJ d A* II 5(A V *) e-^ T ' A *\ (8) 

where ST* is a spanning tree in T*. 

We shall now construct the dual version of this model. We can write 

exp(-AP* 2 (t,)) = ^ J dp{u) exp(-^ + 2ip(t i )P*(t i )). (9) 

Using this integral representation of the plaquette, as well as the standard integral 
of the 5-function: 

r dav e 2ia ' ,A i' 

n <ha*) = / n ^ — (io) 
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the integration over the y4*-fields in T* can be performed, giving rise to a set of delta 
functions: 

n^(a,.+ E P(ti)l (11) 

where the product is over the dual links I*, and «;» = if the dual link I* does not 
belong to the chosen spanning tree. Notice that the sum in the argument of each 5 
function has exactly four terms. Performing the ai* integration in ( |T0D eliminates 
iV 4 (T) - 1 of the N 3 (T) 5-functions in (0). The remaining N 3 (T) - iV 4 (T) + 1 
5-functions can be viewed as the Bianchi identities for the N 2 (T) variables p(ti). 
Since there are N 2 (T) — (N 3 (T) — (A^ 4 (T) — 1)) independent variables p(ti), we can 
implement the Bianchi identities by introducing £7(1) gauge fields A\ living on the 
links of the triangulation T, rewriting p(ti) = P(U) = Independent gauge 

transformations can be done at vertices, and after gauge fixing we are left with 
Ni(T) — (iVo(T) — 1) independent A\ fields, i.e. precisely the same as the number of 
independent p(ti) fields. 

We end up this discussion by writing down the identity valid for each triangula- 
tion: 



/n d A** n ^A,,) e- s ^ T ^ = ff -i^+i^+i e -iE ti inft (12) 

x [UdA l n S(A v )e 

^ 1 ;'e«T 



-S g (T,A) 

I'eST 



where 



S 9 (T,A)=J2^ (13) 

and P(ti) is the oriented plaquette variable (Q) expressed in terms of the fields Ai 
living on the links of the triangulation T. For the partition function (B) we get 



T ^ T •* I I'eST 



-S g (T,A) 



(14) 



where ST is a spanning tree in T and 



3 1 
«4 = — — ln7r, and k 2 = k* 2 — — ln7r. (15) 

Thus, apart from a shift K4 — > k* a — | In 7r and K 2 — > k 2 — \ ln7r, the two models are 
dual up to the weight factor — | hr A- 

The choice of the action parameters A should be made at this point. One can 
argue that it is natural to choose A — ^/ (U), where o(ti) is the order of the triangle. 
In the T* implementation o(ij) is also the length of the boundary of the plaquette 
dual to ti, and, assuming the plaquette to be flat, also proportional to it's area. Since 
P*(ti) signifies the flux of the plaquette dual to t i: the field strength is proportional 
to P*(ti)/o(ti), while the volume element associated with ti likewise is proportional 
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to o(ti). This argument assumes that the plaquette is a flat two-dimensional object 
with a constant field strength. This needs not be the case, so another possible choice 
could be Pi = l/o a (ti) with some power a different from 1, which leaves space for 
some non-trivial geometry of plaquettes. 

The model for a single non-compact Abelian gauge field can readily be general- 
ized to n g copies. For a given triangulation, T, these models will be non-interacting, 
but summing over all triangulations one introduces interactions between the copies, 
mediated by the geometry, as well as interactions between geometry and the matter 
fields. The final action we use can thus be written as 

s >l P ^ (16) 

where the index k numbers the copies of Abelian gauge fields and takes values 
l,...,n g . If a — 1 the model is dual to the original AGM model except for a 
measure factor 

iUEMO (17) 
z t 

for each triangulation and a shift in the gravitational and cosmological constants 
by 

K^-^K^^ ^" m7r an d K 2 K 2 + l n7r - (18) 

We call this version of the model the dual Abelian gauge model (DAGM). 
We can make two remarks at this point: 

• The whole derivation made explicitly use of the fact that the topology of the 
manifold is that of a four-sphere. For different topologies the duality relations 
become slightly more complicated and new degrees of freedom appear, related 
to the possible topologically non-trivial boundary conditions for the gauge 
potentials. 

• Similar discussion can be made in the three and the two-dimensional cases. 
In three dimensions in the T* implementation the dual of a plaquette is a 
link. The duality relates the gauge theory on dual plaquettes to a massless 
bosonic theory, where the bosonic field live on vertices, again up to a ultra- 
local measure term. The T implementation will have it's dual in the form of 
a massless scalar field living in the centers of the simplices. In two dimensions 
the duality transformation permits to integrate out the gauge field completely, 
leaving only the measure term. 



3 The algorithm 

The algorithm describing numerical simulations of the four-dimensional simplicial 
gravity was described in many papers. Here we shall not repeat details concerning 
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the five geometric "moves" which became a standard. There are however some 
peculiarities of the present simulation, which make the updating in the gauge sector 
rather non-trivial. 

In order to perform a Monte Carlo simulation of the model we found it conve- 
nient to work directly with the gauge invariant plaquettes rather than the gauge 
fields themselves. Occasionally we shall need the gauge fields and we shall then 
reconstruct them from the plaquette variables. This is necessary only in order to 
control round-off errors, which eventually will lead to a violation of the Bianchi 
identities. The geometric "moves" create and destroy simplices. In effect also the 
new gauge potentials, located at the interfaces between simplices are created and 
destroyed. Since the gauge action is Gaussian we shall use the heat bath method 
to generate the new variables. In order to find the detailed balance condition it is 
always necessary to compare the "move" and it's inverse. 

The addition of the new gauge potentials does not create any problems and 
can be handled in a way closely resembling that of the Gaussian scalar fields. The 
procedure has two steps: in the first step we decide if the move will be performed. If 
the answer is yes we proceed to the second step, where the new gauge potentials are 
generated from a (multi-dimensional) normal distribution. In the process the old 
plaquette variables are modified and new plaquette variables appear. The Gaussian 
form of the gauge potential guarantees that both the modified and new plaquettes 
will remain restricted in value. The acceptance probability for the first step can 
be expressed in terms of the gauge invariant plaquette variables. In the second 
step the values of new gauge potentials are generated, but only the gauge invariant 
information (the plaquette variables) is stored. 

Situation may be quite different when we try the inverse operation deleting some 
gauge potentials. If we would decide to store the values of the potentials and simply 
subtract these values from plaquettes, the resulting change in the gauge action could 
become arbitrarily large. The reason is that only the gauge invariant combinations 
of the gauge potentials are physical (and in the Monte Carlo simulation restricted by 
the gauge action), while the gauge potentials themselves may become large unless 
some form of the gauge fixing is imposed. Introducing the gauge fixing (a spanning 
tree) is however a global problem and therefore maintaining the gauge condition in 
the numerical simulation, where the geometry is dynamical would be very costly and 
impractical. We have developed a different approach, which is again a modification 
of the heat bath algorithm. In the first step we decide if the inverse move will 
be performed. In order to find the acceptance we integrate over all possible gauge 
choices of the gauge potentials to be deleted. The resulting formula is local and 
expressible in terms of the gauge invariant quantities. If the inverse move is accepted 
we proceed to the second step, where the potentials to be deleted are generated from 
the appropriate normal distribution. After the inverse move only the gauge invariant 
information is stored. 

In the program the situation is slightly more complicated. In order to maintain 
the detailed balance at each move we have to consider the acceptance for both 
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the move and it's inverse. Also the geometric moves used in the four dimensional 
updating of geometry in general are rather complicated and in each move some gauge 
potentials are created, while other are destroyed. We shall report more details about 
the updating procedure in the appendix. 



4 Observables and measurements 

Basic observables in simplicial quantum gravity have already been discussed in a 
number of articles. Since we are here interested in gauge fields it would be natural 
to measure gauge invariant observables. An obvious gauge invariant quantity is 
the correlator Cf^ u (t)F^ p (0)\, where the average is over all geometries with S 4 
topology and with points separated a geodesic distance r. The gauge field itself is 
just Gaussian, so the propagator in flat space is trivially calculated. In quantum 
space-time, where we average over all geometries , it is not clear precisely how the 
propagator will fall off with the geodesic distance. In fact, the motivation for this 
study was to test if the gauge fields had a stronger interaction with geometry than 
the scalar fields investigated so far. In [H] it was shown how to extract the connected 
correlator for matter fields coupled to simplicial quantum gravity. However, in order 
to obtain a non-trivial result of interest for continuum physics it is first necessary to 
identify a point in the combined coupling constant space of gravity and gauge theory 
where we can take an interesting continuum limit. In the following we concentrate 
on this necessary first step. Thus we study the K2, «4 phase diagram of the theory for 
different values of the parameter a in eq. fll6|) . For a given value of Ki there exists a 
critical point k 4 where an infinite volume limit can be obtained. At this critical point 
one can measure geometric quantities like the entropy exponent 7, which describes 
the distribution of baby universes, and one can measure the fractal dimension of 
space-time. As is standard in the measurements of the fractal dimension in simplicial 
quantum gravity studies, one operates with two, in principle independent, definitions 
of the fractal dimension. A "short distance" fractal dimension, dh, measures how 
the volume V(r) of spherical balls of geodesic radius r grows with increasing r as 
long as r is much smaller than the average radius R of the universe: 

V(r) ~ r dh for r < R. (19) 

A "cosmological" fractal dimension is defined by the way "macroscopic" distances 
scale with the cosmological constant or with the average space-time volume. For 
instance 

^)~r*-*F(^), (20) 

where F(x) ~ x b for x — > 0. In the four- dimensional gravity models studied so 
far one has always found that dh = du- Although there exist models [14| of fractal 



structures where dh 7^ dn, one expects that dg = d h if a sensible thermodynamic 



limit exists [12. [T 
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Figure 1: The (k,2, ^4) diagram of the AGM and that of the dual model (DAGM) 
studied in this paper. Error bars are smaller than the point size and the line is only 
to guide the eye. n g = 3 and = 4000. 



Since the transition between the branched polymer phase and the crumpled phase 
of four dimensional gravity is characterized by the appearance of singular vertices 
of very high order, we also study the distribution of the order of vertices, links, and 
triangles in order to characterize the various phases we observe. 

4.1 a = 1 

Since the model (|T6|) for a = 1 is dual to the original AGM model if we include the 
weight factor (|l7]) and make the shift (|l8l), we can as a calibration reproduce the 
K2, diagram of AGM. This is shown in fig. [I] for the number of gauge fields n g = 3. 
We see perfect agreement. As a general statement we also see approximately the 
same values of 7 and dh as observed in the AGM model, thus verifying the results 
obtained in these investigations. It is illustrated in Table 0. Our measurements were 



^2 


7 


d h 


4.0 
4.5 


-0.30(1) 
-0.12(1) 


3.95 (5) 
3.57 (3) 



Table 1: Measurements of 7 and dh for the DAGM model. 



n„ 



3 and N A = 4000. 



made on the N4 = 4000 system with n g 
only dh (and not dn) by using standard techniques 
of 4-simplices at distance r from a given 4-simplex was measured and a fit was made 



3. It was therefore possible to measure 
15 , 14| |. The number riy(r) 
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to the function C(r + a) dh ~ l where a is the so-called "shift", r is the "4-simplex 
distance" which is a measure of the geodesic distance between 4-simplices. The use 
of the shift permits us to estimate dh even for this moderate size of the system as 
can be seen in fig. |2|. A slight modification of the shift a results in a much poorer 
fit for small r. 
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Figure 2: The logarithmic derivative of the two point function ny(r). Here x = 
(r + a)/V 1 / dh where V = For a suitable choice of the shift a a power law 
defining dh can be obtained. 

However, the important point is that for a = 1 and without the measure term one 
moves directly from the crumpled phase and into a branched polymer phase. There 
is no trace of any crinkled phase and the interaction between the gauge fields and 
geometry is as weak as the interaction between scalar fields and geometry observed 
in ||16|| , provided that the gauge fields are coupled to simplicial gravity in a way 
consistent with the way the scalar fields were coupled to the geometry in as 
discussed in detail in the beginning of sec. |2|. 

In table [| we present the results of the computer simulations for the number of 
gauge copies n g = 3 and 6. The value of k 2 is chosen so large that one is well away 
from the crumpled phase. We studied systems with volume iV 4 = 4000, 8000 and 
16000. Our results are summarized in Table 0. 

In fig. [3] and fig. [| we show the computation of du and dh from the scaling 



properties of %(r). du is computed as described in [15], [L4| using the relation 
%(r) = V l ~ l l dn F±((r + a) /V 1 ^ 11 ). We determine du and a from the optimal 
"collapsing" of the %(r) distributions, dh is obtained as for the DAGM model 
described above. We observe that the expected scaling holds very well and that the 
values for dh and dn are consistent with being the same. 

We have also measured the connected part of curvature-curvature and action- 



action correlators (we refer to |L7| for the details). One expects for two observables 
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y 


iV 4 


^2 


7 


dh 


d H 


3 


4000 


4.5 


0.49 (1) 


2.1 (2) 


2.03(7) 




8000 


4.5 


0.51(1) 








16000 


4.5 


0.483(4) 






6 


4000 


7.0 


0.44 (1) 


2.2 (2) 


2.05(6) 




8000 


7.0 


0.47 (2) 








16000 


7.0 


0.470(5) 







Table 2: Measurements of 7, dh and dn when a = l 

A h i = 1,2 that 

G 1 ^(r) = (^)G 11 (r + 5 A J, (21) 
and that the connected part of the correlator is given by 

G^ 2 (r) = G AlA ' 2 (r ) - (Ax)G 1A2 (r + 5 Al ) - (A 2 )G 1Al (r + W (22) 
+ (A 1 )(A 2 )G 11 (r + ^ 1 + ( 5 yl2 ). 

In the above formulas "1" refers to the unity operator and G u (r) = %(r). In 
the case of curvature-curvature correlator we determine the best (Ai) and Si by 
collapsing G 1Ai (r) with G n (r). We use these values in (|2"2]) in order to produce the 
plot in fig. [5]. Similarly we proceed for the action-action correlator. Observe that 
curvature-curvature fluctuations are short range and independent of the system size. 
The same is true for action-action correlators which fall off at the same scale as the 
curvature-curvature ones. Identical plots are obtained for the AGM model so we 
conclude that the measure factor does not smoothen the configurations enough for 
the fields to fluctuate at large scales. The observed behaviour means that using the 
scaled distances the correlations become ultra-local. Similar pattern was observed 
in the crinkled and crumpled phases. This remains in the qualitative agreement 
with the curvature-curvature correlations measured in pure gravity where the 
definition of the connected correlator was different. 

The value for 7 as well as those of dh and dn are consistent with those of ordinary 
branched polymers. The branched polymer interpretation is also corroborated by 
the fact that (N Q ) /N± ~ 1/4 which is the upper kinematic bound expected from 
branched polymers, plus the fact that we see no vertices of high order. 

4.2 a^l 

With the negative results of the above subsection in mind we turn to the more 
general action (|16"D with a 7^ 1 in a search for new and potentially interesting fixed 
points from the point of view of continuum physics. By performing the duality 
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Figure 3: Calculation of du using finite size scaling of riy(r). The left figure cor- 
responds to n g = 3, k 2 = 4.5, d H = 2.03 and the right one to n g = 6, k 2 = 7.0, 
d H = 2.05. 

transformation for a general a we obtain an action 

£(2>W(M))+^lno(t). (23) 
t k=i / 

Thus there is a possibility that a choice of negative a might bring us to a phase 
resembling or being identical to the crinkled phase, since the negative weight factor 
will imitate the situation encountered in the AGM model. We show the result of 
the measurements for the number of gauge fields n g = 3 and a = —0.5 and —1. The 
value of k,2 is chosen such that we avoid the crumpled phase. 

The results of the measurements of 7 and dh is shown in Table [| There is a clear 
tendency for the values to 7 and dh to drift from the branched polymer values toward 
the crinkled values, in accordance with expectations. Some limited statistics that 
we obtained for the a = —1.5 model suggest that 7 becomes even more negative. 
We have taken this as evidence that we see not really new physics with the modified 
weights, but rather an effective change in measure, more or less in the same way as 
the original AGM model differs from our model by measure term. 

5 Discussion 

Gauge fields contribute far more to the conformal anomaly than scalar fields. Thus 
they are expected to play an important role in effective models of four-dimensional 
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Figure 4: Calculation of dh for n g = 3, k 2 = 4.5 (left) and n g = 6, tz 2 — 7.0 (right). 
Here we use x = (r + a)/V 1 ' dh . The shifts a are —0.5 and —1.5 respectively. 



a 


7 


dh 


-0.5 
-1.0 


0.30(2) 
-0.51(1) 


2.65 (3) 
3.5 (2) 



Table 3: Measurements of 7 and c4 for a < 0. n 9 = 3, A/4 = 4000 and k 2 = 4.5. 

quantum gravity where the (infrared) dynamics is dictated by the conformal anomaly 
pj. The conjectured scaling behavior of these model are quite similar to the observed 
(pseudo) scaling behavior observed in four- dimensional simplicial quantum gravity, 
and this led to the expectation that gauge fields might couple stronger to geometry 
than the scalar fields used so far. So far, computer simulations have revealed only a 
weak coupling between the scalar fields and the geometry, and the presence of scalar 
fields have not led to any quantitative change in the phase diagram of simplicial 
quantum gravity. 

The first Monte Carlo simulations of simplicial quantum gravity coupled to gauge 
fields seemed to be in accordance with the above philosophy. Indeed, for the first 
time one observed a genuine back-reaction of the matter fields on the quantum 
geometry. However, as already remarked in the original article and verified in 
detail later 0, major part of the interaction between geometry and the matter could 
be accounted for by the term 

Se ff = f EHo(t)), (24) 
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Figure 5: The connected curvature-curvature (left) and action-action (right) corre- 
lation functions. n g = 6 and k 2 = 7.0. The n g = 3 plots are falling on top of the 



n g = 6. 



where o(t) denotes the order of the triangle t. This term is an ultra-local measure 
term and it is (as already remarked in ||) unlikely that it contributes to the infrared 
dynamics related to the conformal anomaly in the model proposed in [0, ^ . However, 
it still left us with the puzzle why gauge fields seemingly interact so much stronger 
with geometry than ordinary scalar fields. 

In the work presented here we have verified the original results obtained in the 
AGM model by an independent Monte Carlo simulation working with gauge fields 
on the lattice dual to that of the of the AGM model. In addition we have shown that 
there is no major difference between the interaction of scalar fields with geometry 
and the interaction of gauge fields with geometry. In fact, if the gauge fields are 
coupled to geometry in a way consistent with the way the scalar fields were coupled to 
geometry, the gauge and the scalar fields have similar weak coupling to the geometry. 
We have explained how the weight factor (|24]) arose in the transformation from gauge 
fields living on the dual links of a triangulation to gauge fields living on the links 
themselves. 

Turning the arguments around, our results lead to the prediction that scalar 
fields living on the vertices of four-dimensional triangulations will couple as strongly 
to the geometry as the gauge fields in the AGM model. 

We extended that gauge field model coupled to geometry by modifying the weight 
given to each plaquette, in order to explore further the phase-space structure of 
matter coupled to geometry. However, the only new phase we could identify is 
presumably identical to the crinkled phase already found in the AGM model. The 
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weight generated by the transformation from the dual links to the ordinary links of 
a triangulation gave support to this interpretation. 



Appendix 

All moves require modification of some plaquettes by introducing new dual links 
into plaquettes or by destroying the dual links. Finding probabilities for the moves 
resembles the problem of the Gaussian scalar fields, however the gauge invariance 
of the action requires some important changes to maintain the detailed balance and 
at the same time to keep the acceptance reasonable. In this section we shall discuss 
the case of one gauge potential and the parameter a — 1. 

Let us first describe the process of introducing new degrees of freedom. For 
simplicity let us discuss only the situation of the move, where a new vertex is added 
in the center of a simplex. There are now 10 new gauge potentials. We consider only 
the part of the gauge action involved in the move. Before the move is attempted we 
have 

10 

exp(~S A (pi)) = exp(- PiPi), ( 25 ) 

i 

where $ = l/o(tj) and pi are the external plaquette variables. The modified action 
becomes 

10 

exp(-S B (p u Xi)) = exp(-J2Pi(Pi + x i) 2 -P'J2 x i M ij x j) ( 26 ) 

i i,j 

where Xj are new fields, f3[ are the new weights, typically j3[ = l/(o(tj) + 1), the 
matrix contains the geometric information about ten new internal plaquettes 
formed in the move. These new plaquettes IT^ can be written as 

n fc = E UkjXj (27) 

M^ = Y,U kl U k] (28) 

k 

where Uki are ±1 or zero reflecting the orientation of the Xi fields. The internal 
plaquettes have always the order three, so (3' = 1/3 typically. When performing the 
move we shall have to 

• decide whether we make the move or not and 



• if yes - generate the new fields with the distribution resulting from (|2 
Let us define 

N B (jPi) = jY[dxiexjp(-SB(Pi,Xi)) = 



7T 10 



\ det(Q) 



ex P (-S B ( Pt )), (29) 
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where Sb{Pi) is the global minimum of Ss{pi, x%) with respect to Xi. This minimum 
can of course be explicitly calculated in terms of Pi and the matrix Q -1 : 

SB(ft) = £#P?-Eft<#Pi- (30) 

i ij 

The quantity wsipijXi) = exp(— Ssipi, %i))/NB(Pi) is a normalized Gaussian 
distribution of the Xi fields. Let Vb = t^b-^b be the probability that the move will 
be performed and, if accepted, the new fields be generated with this distribution. It 
is easy to check that the transition probability P(A — ► B) satisfies 

exp(-S A )P(A -> B) = exp(-S A )V B w B (Pi,Xi) (31) 

= 7r B exp(-S A - S B )- 

The form ([H|) will be useful in order to discuss the detailed balance condition. 
The important point is that the only non-trivial quantity, Ns{Pi) can be calculated 
also in the new configuration {Pi = pi + £j,IT fc }. Using ( |2"§D we can express Sb as 
a global minimum of J2ifli(Pi + Hi) 2 + P'J2k(J^-k + HjUkjVj) 2 with respect to the 
shifted variables y^. The only difference is the shift of variables and the value of the 
minimum is unchanged. This property is very important, because we shall have to 
calculate Mb both for the move and it's inverse. 

Let us now discuss the inverse move: we destroy 10 internal plaquettes 11^ and 
modify values of 10 external plaquettes, eliminating 10 gauge fields produced before. 
The naive proposition would be to store the gauge fields x« and simply to delete them. 
The gauge fields however have no gauge invariant values so the transition performed 
this way would depend on the particular gauge choice - in fact even if plaquette 
variables are small (because of the Gaussian weights) the gauge fields may become 
large and in effect the transition may be blocked completely. 

In the program we use a different approach. Rather than storing the gauge 
fields we try to reconstruct them when needed. When trying to delete fields Xi we 
have first to decide what gauge choice we make. In other words we shall study the 
change Pj — > pi = Pi — Xi where x« satisfy the set of gauge-invariant constraints 
UkiXi = Ilfc. Notice that not all these constraints are independent because of 
the Bianchi identities between the ten plaquettes Ufc. Only six of them have to be 
used, reducing the number of independent degrees of freedom to four. As before let 
us consider only the part of the action engaged in the move. Before the move is 
attempted we have 

10 10 

exp(-5 B (P l ,n fc )) = exp(-£# P 2 -13' £ n fc) (32) 

i k 

After the move we have: 

10 

exp(-S A (Pi, x^) = exp(- £ - x^ 2 ), (33) 

i 

As before we shall have to 
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• decide if we perform the move and 

• if yes, perform it, i.e. choose some values of Xi following from (33|). 
Let me define as before: 

.10 6 

J\f A (Pi,U k ) = I Y[dxi Y[S{Il k - ^2U kj x j )e^,{~S A {P i) x i ) 1 



k 



I 7T 4 

detN 



ex V (-S A (P h U k )), (34) 



where the matrix N depends only on $ and S A (Pi, life) is the conditional minimum 
of S A (Pi,Xi) with gauge-invariant set of conditions imposed on the Xi fields. The 
integration in (0) can be viewed as integration over all possible gauge choices for 
the fields Xi. It is simple to give the explicit formula for S A in terms of Pi, U k and 
the matrix N _1 . We shall not write it here. Let us however note that J\f A (Pi,U k ) 
can be also calculated in the new configuration {Pj,IT fc } — > {pi,0} using a simple 
relation M A {P h U k ) = Af A (p t , 0). 

The quantity w A (Pi, U k , Xi) = exp(—S A (Pi, Xi)/J\f A (Pi, H k ) is a normalized prob- 
ability for Xi satisfying the gauge-invariant constraints. As before let us choose the 
probability V A = n A J\f A (Pi,U k ) to be the probability of performing the move. If 
accepted we generate four independent Xi from the distribution w A . The remain- 
ing six can be calculated using the constraints. We can check that the transition 
probability P(B — > A) satisfies 

exp(-S B )P(B -> A) = 7r A exp(-S A - S B ). (35) 

Equations (|3l|) and (^) can be used to write the detailed balance condition: 

exp(-S A )P(A -»■ B) = exp(-S B )P(B -> A) 

meaning 

^ = 1 

k a 

or v2 = W A (36) 

As usual, this form of the detailed balance condition suggests that in order to max- 
imize the acceptance we need to know both Mb and A/a -The generalization to the 
case of more gauge fields and a ^ 1 is trivial. 

The same method can be used to calculate the weights for other moves. The 
general rule is that we need to consider at the same time both the move and it's 
inverse. 
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